Non-exotic theory of l/f noise as a trace of 
infralow-frequency fluctuations 

A. Ya. Shul'man 

Institute of Radio Engineering and Electronics, Russian Academy of Sciences, 
Moscow 125009, Russia. E-mail: ash@cplire.ru 

Abstract. This report is aimed at reviving the explanation of flicker-noise observations as the result 
of spectral measurement of very low-frequency but stationary narrow-band fluctuations named 
as infralow-frequency noise (ILF noise) [A. Ya. Shul'man. Sov. Tech. Phys. Lett. 7, 337 (1981), 
Sov. Phys. JETP 54, 420 (1981)]. Such a kind of the spectrum analyzer output takes place when 
the ILF-noise correlation time is much longer than the analyzer reciprocal bandwidth. This result 
is valid for both analog and digital spectral measurements. The measured signal is proportional 
in this case to the mean square and not to the spectral density of the noise. The equilibrium 
temperature fluctuations and the defect-motion as mechanisms of l/f noise in metal films are 
reconsidered from this point of view. It is shown that the ILF-noise approach allows to remove 
the main objection against temperature fluctuation model and difficulties within the defect-motion 
model widely discussed in literature. 
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1. INTRODUCTION 

It was shown in Jl|]-||2] that flicker-noise observations can be explained as the result 
of spectral measurement of very low-frequency but stationary narrow-band fluctuations 
named as infralow-frequency noise (ILF noise). Such a kind of the spectrum analyzer 
output takes place when the measured noise spectrum concentrates on frequencies far 
lower than the operation range of the spectrum analyzer and the ILF-noise correlation 
time is much longer than the analyzer reciprocal bandwidth. This result is valid for both 
analog and digital spectral measurements. The measured signal is proportional in this 
case to the mean square rather than the spectral density of the noise, and does not reflect 
in any way the dependence of the latter on the frequency. Within the scope of this model 
the universality of the l/f noise follows from the universality of inevitable fluctuations 
of intensive thermodynamic variables like the temperature and a power-low falling of 
narrow-band filter spectral distributions in the wings. It is of importance here the more 
isolated is the sample under study from the environment the more low-frequency are 
corresponding fluctuations. Thus, it is impossible to reach the low-frequency saturation 
region of the noise spectral density if the duration of the measurements is less than 
response time of the sample - environment coupling. 

The necessity to find not only slow fluctuation process but also with l/f -spectrum 
serves as the reason to call in question some physical mechanisms that by their qualita- 
tive attributes and quantitative estimations could be a cause of the excess current noise. 
Example of this is the temperature model of l/f noise in metal films suggested by Voss 



and Clarke on the basis of convincing experiments |0]. Taken initially with hopes and 
an enthusiasm [4]-[5] this theory has later been rejected as an explanation of 1/ f noise 
in metal films [6],|7]. In addition to the proof that the heat conductivity equation does 
not lead under no circumstances to the 1/f -spectrum of temperature fluctuations fSD, 
another weighty argument against this model was the absence of spatial correlations of 
the current fluctuations [9], [10]. 

On the contrary, there is no need in the scope of the ILF-theory for exponentially- 
wide distribution of relaxation times of any kind |0]. The interpretation of the spatial 
correlation measurements is also to be reconsidered since the mean square of the conter- 
minous fluctuations is measured in the case of ILF-noise rather than the spectral density 
of them. It will be shown here that all abovementioned objections against temperature 
fluctuations are withdrawn from the ILF-noise point of view. 



2. THE NOISE SPECTROSCOPY AND INFRALOW-FREQUENCY 

FLUCTUATIONS 

2.1. Analogue measurements 

It was shown in lH-H that the output G out (a>o) of an analogue spectrum analyzer as 
a function of the resonant frequency (Do can be described in a good approximation by 
the convolution integral 

oo 

G m,t (coo) = I J d -^X (^^j G y (a), A/ao « 1, < 1^1 < <*>, (1) 

— oo 

(the notation is slightly changed). Here G y ((0) is the spectral density of the stationary 

random input process y(t) with (y 2 ) < <*>, ( C0 ° A C ° ) plays the role of the instrumental 
function of the spectrum analyzer normalized by the condition J?T(0) = 1, A is the 
effective bandwidth of the analyzer defined by 

oo 

CO 

An important point is that Jt{(o) and G y ((0) are nonnegative and normalizable 
functions of the frequency and can be approximated by 8 -function each at proper 
interrelation between the effective bandwidths. If J^(co) is noticeably more narrow- 
band function then 

G out ((O ) = G y (co ) 

and the result of measurement G out (coo) is the common- adopted estimation of the true 
spectral density G y (oo). On the other hand if the investigated spectrum has a narrower 
contour than the instrumental function, then the result of the measurements will be the 
contour of the instrumental function 

G-(^) = ijr(^^)<v 2 >, (3) 



where co m is the frequency at which the maximum of G y (co) takes place. Spectra are 
always measured with a narrow-band filter, when the quality factor or the figure of merit 

Q = /o/A»l,(/ = fflb/2tf) (4) 

. It is now obvious from Eq. © that when the spectral density of the noise is narrower 
than J^(ft>), and the noise frequency is low enough co m <C C0q [this is the definition of 
the term "infralow-frequency" noise], then in the usual Q = const (coq) measurements 
the output should have a dependence of the l/f type: 

G m "((Oo) = Qj*r(27iQ) (y 2 ). (5) 
Jo 

On the other hand, if the measurements are carried at A = const (coq), then the output 
signal has the frequency dependence of the far wing of J^(co), which likewise has as a 
rule the form 1/ f n . In either of these case the measured values are proportional to the 
mean square (j 2 ) of the input random process. 



2.2. Digital spectral analysis 

In some respect the analysis of digital methods, from our point of view, is even 
simpler, since the formula similar to Eq. ([TJ is already known. There is the exact 
expression of the form 

CO 

G out (coo) = T J ^W(T o ((0o-(o))G y ((o), < coq < °°, (6) 

— CO 

where Tq is the time slice of the one realization of the random process and W(co) is the 
spectral window. Widely used windows can be written in the form 



W{(0) 



sin ( (o To/n) 
(oTo/n 



(7) 



where n = 2,4 corresponds to the Bartlett and Parzen window, respectively. 

The account for finiteness of the time step To and application of discrete Fourier 
transformation lead to a little bit changed expressions for spectral windows. For Bartlett 
window one has ||2j] 

Af 2 snr(ct)To/2) 

and Tq = Ntq, N is the number of data points in the measured realization. So long as 
the frequency co is small compared with the Nyquist frequency, defined by the condition 
(On = 2nfy = tz/xq the continuous and discrete Fourier transformations yield practically 
identical expressions for the instrumental function (spectral window). However, since 
Wd is the periodical function of CO, the formula © can give rise to the aliasing effect. 
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FIGURE 1. Dependence of the spectrum analyzer output on the tuning frequency, Eq. (11), under 
condition of the constant quality factor Q = /o/A. Two typical values of Q are presented corresponding 
to the equivalent noise bandwidth of 20% and 5%, respectively. The spectral density of the noise at the 
analyzer input, Eq. (10), is also shown for comparison. 



To avoid this, the input process G y (co) should be filtered by low-pass filter with cut-off 
frequency of order of fy. With this condition in mind, the comparison of Eqs.lH])-© to 
Eqs.©-© shows that the effective bandwidth is A = I /To in the digital spectral analysis 
with the Barttlet window. 

Comparing Eqs.©-© with Eqs.Q and © we can easily verify that in the case of 
digital processing of ILF noise measurements, an erroneous observation of the spectrum 
of the 1/ f n type will also take place. The digital measurements are always carried out 
with a constant bandwidth (~ l/7b), at least in each of several frequency intervals into 
which the entire investigated band is subdivided. As a result of effective "averaging" of 
the resultant curve over all the sub-bands, a spectrum of the type 1 //" can be obtained 
with 1 < n < 2 for Barttlet window under conditions of ILF-noise measurements. 



The observable signature of ILF noise measurements is shown in the figure [T] It 
was calculated in the case of Lorentzian form of both the noise spectral density and 
the narrow-band input filter of the spectrum analyzer. The following expressions were 
substituted in Eq. CQ) as the J^(co), A and G y (co) 



2.3. An illustrative example 



f/4 




b 2 - f/4, A = y/4 



(9) 



(a>i - (o) 2 + y 2 /4 




(10) 



to obtain the explicit formula for the spectrum analyzer output 



ft> 1 2 + (v + 7/2) i 

The parameters C0q and 7 are defined by the equation of the simple resonant filter 1 

x + Jx + COqX = y(t). 

It is easily seen that the condition © for the analyzer to be narrow-band (in our case it 
takes the form ny/lcoo <C 1) still does not determine uniquely the G out ((0\) dependence. 
It is also important to know the relation between 7 and v, i.e., between the widths of 
the input spectrum and the resonant filter. At y/v <C 1, i.e., in the case of a relatively 
narrower filter, we have the usual result G out ((Q\) = G y ((Q\). In the other limiting case, 
however, when 7/ v ^> 1 , we obtain 

G-(«0 = <^ 3 ^>g (12) 

From Eq. (fT2l it follows directly that if the measurements are carried out in this 
situation at a constant Q, i.e., under the condition y/(0\ = const(coi), then the output 
signal will imitate the presence of noise with a spectrum 1 // at the input of the analyzer. 
It is clearly seen in Figdthat although at high frequencies (ft>i ^> v) the true spectral 
density of the noise (ITOb is everywhere proportional to 1 /ft) 2 and the condition © is 
satisfied, all this is still not enough to let the region of the maximum Jf{(0\ — ft)) to 
make the main contribution to the integral in Eq. ([T]). Conversely, as soon as the filter 
bandwidth 7 becomes larger than the width of the input spectrum v, a transition from 
the expression (flOb to (fT2l takes place quite rapidly. In other words, the output signal 
of the analyzer ceases to yield information on the form of the spectrum of the input and 
begins to vary with frequency in accordance with the 1 // law. 

The measured curves of such a kind can be found in literature without a discussion of 
possible ILF-noise interpretation. As a result, some speculations take place in order to 
explain the observed transition of measured spectral density from 1/ f 2 to 1/ f either by 
means of an unusual frequency dependence of the sample properties (see, for instance 
lUTlo or a hypothesis is offered based on a decomposition of measured noise spectrum 
on the sum of Lorentzian and 1/ / components originated by different causes lfl2ll . 



3. 1 JF NOISE IN METAL FILMS 

Investigations of the excess current noise in metal films purposed to shift the study of 
unclear l/f noise onto the objects that should be physically more defined than non- 
uniform semiconductors or thermionic emission. If refer a start of that researches into the 



To compare Eqs. (11) and (12) with corresponding Eqs. (10) and (12) in [2] one has to correct the 
misprint there and change 7 — » y/2 taking into account the condition (4). 



beginning of 1970th, then it is to be ascertained that more 30 elapsed years and numerous 
new obtained data have not led to the decision of a problem of origin of such a spectrum. 
Although observed dependencies of spectral density of resistance fluctuations at any 
fixed frequency on various parameters of a sample material or measurement conditions 
point out as if these variables are related with 1/ / noise, it is failed to deduce from here 
such its statistical properties as the spectral form or the behavior of spatial correlations. 
The approach from the point of view of ILF-noise allows to remove an acuteness of 
some contradictions and restate the problem for additional experimental studies. 



3.1. The temperature fluctuations 

Temperature fluctuations as the cause of 1 //-noise have been rejected by three prin- 
cipal arguments |0]-|0]: 1. The impossibility to deduce the 1 // spectral density for fluc- 
tuations of the volume- averaged film temperature. 2. The absence of correlations of 
current fluctuations in spatially dividual parts of the film. 3. An effect of the substrate on 
the measured resistance fluctuations. Let us consider these conclusions from ILF-noise 
point of view. 

To analyze the Voss and Clarke experiments it is necessary to address to Eqs.© and 
®. Assuming G y ((6) is ILF-noise and co is not too close to Nyquist frequency (On we 
have for the relative voltage fluctuations (2D 

G^jf) _ 2sin 2 (^T /) P 2 (8T 2 ) _ 1 /3 2 7/ 2 

£/2 n 2 (f/A) f n*(f/A)3N a f { ) 

Here G + (f) = 2G(f) is the spectral density in / > frequency semi-axis, /3 = 
dlnR/dT, R is the resistance of the film, (8T 2 ^ = ksT 2 /C is the mean square of the 
equilibrium temperature fluctuations averaged over the film volume, C is the film heat 
capacity. In the case of the room temperatures and monatomic metals C = 2>kgN a and N a 
is the number of atoms in the film. The substitution of experimental data [3] combined 
with the hand-book material parameters for the Bi and Au films in Eq. (fT3l and hand- 
made estimations of Ag ; = 2Hz and A^ u = I Hz results in[2] 

G^ ut (10Hz) = 15 .10~ 17 Hz _1 , G^° ut (10 Hz) = 0.58 -10- 17 Hz~ l 

The measured values are 13 • 10 l6 Hz 1 (Bi) and 0.6 ■ 10 16 #z 1 (Au). The differ- 
ence just in 10 times for both cases suggests that there might be some problem with 
conversion of digital Fourier-transform coefficients in 1 /Hz spectral density. Anyway, 
the obtained of-order-of-magnitude coincidence is worthy of notice, taking into account 
the lack of fitting parameters and indeterminacy of the film dimensions mentioned in 
10]. The 1/ ' f a dependence of the analyzer output with 1 < a < 2 follows from Eq. (fT3l) 
as it is explained in Sec. 2. 2. 

As for the lack of spatial correlations of the volume- averaged variables, it is necessary 
to note, first of all, that the dependence of the relative correlation function like in Eq. (fT3l) 
on the sample volume CI in the form of 1 /CI may occur only if the characteristic scale r c 
of the spatial correlation function (8T(r,t)8T(r\,t\)) is much less than the dimensions 



of the sample. The discussion of this point can be found in El from general point of 
view and in [Jg] with reference to 1 //-noise problems. Therefore, the spectral density 
of the frequency Fourier components should be spatially correlated if they are formed 
from the spatial Fourier components with characteristic wavelength of order of the film 
dimensions. Thus there is no sense to look for the spatial correlations between different 
parts of the film if the I/O. dependence of correlation function is observed and vice 
versa. 

This receipt regards to the case of the usual spectral measurements. In the case of 
ILF-noise the output signal is proportional to the mean square of temperature fluctua- 
tions taken in the same instant time. In an equilibrium macroscopic body out of phase 
transition region such fluctuations are always uncorrelated on macroscopic scales in ac- 
cordance with the principle of additivity of extensive variables in thermodynamics. The 
explicit prove of this statement for the spatial-averaged temperature fluctuations can be 
found in [8]. So in the case of ILF-noise the absence of spatial correlations can not be 
an argument against of the temperature-fluctuation origin of l/f -noise. 

To consider the dependence of l/f noise on the substrate supporting the film and on 
a degree of the film-substrate thermal contact one need to account for the last stage of 
the time relaxation of the averaged temperature fluctuations, i.e. a heat interchange with 
the environment. The Newton boundary condition for the thermoconductivity equation 
gives rise to the temperature relaxation with one characteristic time T, resulting in the 
simple expression for the spectral density of the form 

G T ((D) = (8T 2 ) - 2T 7 . (14) 

As T is very large in comparison with the reciprocal lowest measured frequency the 
ILF noise should take place and the dependence on the substrate in Eq. (TT31) comes from 
the heat capacitance C, which is determined by all structure under study including an 
involved part of the substrate. The latter depends also on the heat-transfer factor of the 
film- substrate interface and can increase as X decreases, giving rise to the decrease of 
G out (co). As 1/t gets into the measured frequency range the usual spectral measure- 
ments takes effect and l/f behavior of the spectral density can give place to an ordinary 
one. 

So both characteristic dependencies of spectral measurements on the substrate, which 
were discussed in H without definite conclusions, can be understood on the base of 
ILF-noise model within the scope of the temperature fluctuations' theory as well as two 
other problems considered in this section. One can believe it is premature to discard the 
temperature fluctuations as the source of l/f noise in metal films. 



3.2. The defect-motion model 

After Eberhard and Horn have found out the activation dependence of the l/f -noise 
on the metal film temperature ifLHl the opportunity of 1 / /-noise generation by defects 
of the lattice was taken into consideration. As usual in the field of flicker noise there is a 
contradiction between the observable well-defined activation energy and the necessity 



to have a wide-spread range of these energies to produce l/f spectral dependence 
over an extended frequency range. It was note in 111411 there is no possibility to avoid 
this contradiction using only one activation process. It was suggested to divide the 
defect contribution into the current noise by two parts connected with the fluctuations 
of the defect number due to a vacancy creation/annihilation and the fluctuation of the 
free-carrier- scattering cross-section due to a defect motion with different activation 
energy. The last process should be giving the bread distribution of the activation energy. 
However, the observed noise activation energy in various metal films are less than known 
vacancy formation energy approximately to 10 times 11411 . Hence the observed activation 
energies should be related with the motion of defects, leaving no place for the usual 
construction of 1 //-noise model. 

Kogan and Nagaev [15] suggested to connect the motion-defect contribution to the 
flicker noise with some characteristics of the internal friction process. The reasonable 
estimation of possible noise level was obtained for that mechanism (see also $7$) but the 
origin of the wide distribution of the activation energy remained unexplained [6]. The 
ILF-noise model does not required any real l/f spectrum of fluctuations. It should be 
sufficiently low-frequency only. The internal friction processes rather meet this require- 
ment having characteristic relaxation frequencies in the diapasons of 10 12 , 10~ 8 , 10 
Hz with magnitudes of reciprocal quality Q m l x in the 10 _1 -j- 10~ 2 range L 1 6] . One has 
to recall that the value of Q ~ 10~ 4 for the background level of the internal friction 
was used in [15] with a positive conclusion. 



4. CONCLUDING REMARKS 

It is worthy to mention the work lfl7ll where the magnitude of l/f noise in SQUIDs 
was evaluated as a result of ILF- fluctuations of the current in superconducting or normal 
metal shield surrounding the SQUID. The results obtained there keep their significance 
for high-temperature SQUIDs, too lITsll. 

It seems evidently that the ILF-noise model could remove contradictions in existing 
theories of l/f noise and should be involved in the analysis of experimental results in the 
field of flicker-noise investigations. A necessary change in the experimental technique 
has been shortly discussed in yj]. 
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